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Abstract 

A simple radiation condition at infinity for time-harmonic massive 
Dirac spinors is proposed. This condition allows an analogue of the 
Cauchy integral formula in unbounded domains for null-solutions of the 
Dirac equation to be proved. The result is obtained with the aid of meth- 
ods of quaternionic analysis. 
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1 Introduction 

There are two principal ways to guarantee the unicity of solution of a not homo- 
geneous differential equation of mathematical physics in the whole space. The 
first is the famous limiting absorption principle which in the case of some oper- 
ators, for example the Hclmholtz operator, has a very natural physical meaning 
reflected in its name (see, e.g., p3). The medium is assumed to be slightly 
absorbent with a corresponding absorption parameter e. For e > the inho- 
mogeneous equation (under some additional natural conditions) has a unique 
solution u e decreasing at infinity, and the solution for e = is obtained as the 
limit of u e when e — > 0. This approach, based on a deep understanding of the 
physical nature of the differential operator, is not applicable in all cases. An- 
other approach is to impose a radiation condition at infinity which "prohibits" 
the existence of some unwanted solutions. Such a prohibition is related to the 
unicity of an integral representation for null-solutions of the differential opera- 
tor. In the case of the Hclmholtz operator the corresponding radiation condition 
was introduced by A. Sommerfeld [|10|. 
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For the Maxwell equations such condition was proposed by C. Muller S and 
S. Silver ||. Recently, A. Mcintosh and M. Mitrea in JtJ introduced a radiation 
condition for a special perturbed Clifford-Dirac operator. We use their approach 
to obtain a radiation condition for the operator D+M" , where D is the Moisil- 
Theodoresco operator and M is the operator of multiplication by a complex 
quaternion ~ct from the right-hand side (for exact definitions see Section 2). 
Then using a simple matrix transformation proposed in which converts the 
operator D + M into the Dirac operator in its traditional form through the 
7-matrices, we rewrite the radiation condition for the Dirac spinors in a form 
suitable for applications and obtain an integral representation for null-solutions 
of the classical Dirac operator in unbounded domains, based on the proposed 
radiation condition. 



2 Preliminaries 

2.1 Quaternions 

A quaternion a can be represented as a linear combination of its components 
({a/c} C R for real quaternions and {etfc} C C for complex quaternions) and the 
elements of the orthonormal basis ik as follows 

3 

a = 2J akik 

k=Q 

where io :— 1 and {ik \ k = 1, 2, 3} are the quaternionic imaginary units, which 
have the following properties: 

if, = —1, k = 1,2, 3; 



«1«2 



"«2«1 = »3; «2»3 



"«3»2 = ll] 



«3«1 



-«1«3 = ^2■ 



(1) 



The imaginary unit in C is denoted as usual by i and by definition commutes 
with quaternionic imaginary units. 



1 ■ ik = ik ■ J, 



0,3. 



In some cases it is useful to represent a quaternion a as the sum of a scalar 
part Sc(a) = ao and a vectorial part Vec(a) = ~ct = a\i\ + a^ii + a^i^. Using 
this representation we can define the product of two arbitrary quaternions a 
and b in the form 



a ■ b = aobo 



b ) + a b + ~a*b + 



~a? x b 



(2) 



where (•, •) is the scalar product and [• x •] is the vector product. 
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The quaternionic conjugation of a quaternion a = oq + ~ct is defined as 
a — ciq — ~~ct . From (|^) it is easy to observe that 

a ■ a = <2g + a\ + a\ + a 2 =: \a\ 2 , (3) 

and to obtain an important property of quaternionic conjugation: 

a ■ b = b ■ a 

Note that \a ■ b\ = \a\ ■ \b\ . 

One important conclusion which can be made using ([}]) is that each non-zero 
real quaternion a is invertible and its inverse is given by 




We denote the set of zero divisors in H(C) by 6: 

S := {a e H(C) \a^0;3b^0:a-b = 0}. 



Lemma 1 ([^]) (Structure of the set of zero divisors) Let a 6 H(C) and a ^ 0. 
The following statements are equivalent: 

1. ae&. 

2. a-a = 0. 

3. al = "a" 2 . 

4-. a 2 — 2a^a = lit a. 

As the modulus introduced by (||) does not in general give information about 
the absolute values of the components of a complex quaternion when it is a zero 
divisor, another kind of modulus is frequently used 




where |afe| 2 = a^a^ and "*" stands for the usual complex conjugation. (Q) can 
be expressed also in the following ways 

| a. | ^ = |Rc a\ 2 + |Im a\ 2 = Sc(a • a*) = Sc(a* • a). 
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2.2 Quaternionic differential operators 



Let CI C R 3 and / E C x (fi; H(C)). The Moisil-Theodoresco operator is defined 



as 



k=l 

where dk '■— An important property of I? is obtained from (^|) 

D 2 = -A, (5) 
where A = df + 9| + d 2 is the usual Laplace operator. 
We will consider the following operator 

D± v :=D±vI 

where v is a complex number, which was studied for the first time in [|| (for v 
real) . With the aid of the operator D v it is possible to factorize the Helmholtz 
operator 

A + v 2 =-(D + v)(D-v) =-D v D- v . (6) 
Suppose that Q v is a fundamental solution of the Helmholtz operator 

(A + v 2 )Q v = S. 
Then, using (JsJ) we obtain the functions 

K± v :=-{Dtv)Q v (7) 
which are the fundamental solutions of D± Ul that is 

D± V K± V = 5. 

Let us assume that Imv > 0. Under this condition, the fundamental solution 
of the Helmholtz operator is chosen as 

p iv\x\ 

e v ( x ) = -—-, (8) 

47T \X\ 

in which case it fulfills the Sommerfeld radiation condition at infinity and has 
a clear physical sense of an outgoing wave generated by a point source at the 
origin. 

Substituting (ph into equality (pi), we obtain 

/ x \ 

K.± v (x) = —D6 v {x) ± vd v {x) — ±v H g ~ i v T~\ 9v( x )i 

\ \x\ F / 
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where x ■=Ylk=i x k^k- These functions were obtained in (see also |6[ Section 
3]). 

Let us assume that f2 is a bounded domain in R 3 with a piecewise smooth 
boundary T := dfl. 

Theorem 2 ( Quaternionic Stokes 'formula) Let f and g belong to C 1 (f2; H(C))D 
C(0;H(C)). Then 

[ ((D r f(y))g(y) + f(y)(Dg(y)))dy = [ f{y)jt{y)g{y)dT yi (9) 
Jn Jr 

where ~rt(y) is the quaternionic representation of the outward unit normal to 
the surface T : Tt{y) = J2k=i *fc n fc(y) an d D r f '■= J2k=i®kfik is the right 
Moisil-Theodoresco operator. 

The proof of this fact can be found in || Chapter 4] or ]|, p. 86]. 
Let us consider the following integral operator 

K v [f]{x):=- j K u {x-y)!t{y)f{y)dT v , x e R 3 \ T (10) 

which can be considered as an analogue of the Cauchy integral operator from 
complex analysis due to the following fact. 

Theorem 3 (Quaternionic Cauchy integral formula) Let f 6 C 1 (fl; H(C)) n 
C(Q;H(C)) and f e ker D„(fl). Then 

f(x) = K v [f](x), VxeO. (11) 



2.3 Dirac operator in quaternionic form 

Denote $ := Q(t, xi, X2, —X3). The domain Q is assumed to be obtained from the 
domain Q C K 4 by the reflection xz — * —2^3. Let us consider the transformation 
denoted by A and introduced in the following way. A function $ : Q C R 4 — > C 4 
is transformed into a function p : Q c R 4 — > H(C) by the rule 



p = := - (^-($1 - <J> 2 )«o + «($o ~ *3)«i - ($0 + $3)22 + + $2)«3 
Then the definition of the inverse transformation A^ 1 is 

A^[p\ = (-ipi - fc,-po - ip3,Po ~ ipz,ipi ~ fa)- 
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These transformations can be represented in a more explicit matrix form: 



p = m = 

and 





/ 





-1 


1 






1 $ \ 


1 




i 








—i 




Si 


2 




-1 








-1 








V 





i 


i 


o / 




V *3 / 




f 





—i 


-1 


N 




( * \ 






-1 








—i 




Pi 






1 








—i 




P2 









i 


-1 


) 




V ?3 / 



The Dirac 7-matrices in the standard (Q, Dirac-Pauli form are 



7o := 



72 := 



/ 1 


\ 








/ 








-1 


\ 





1 











- 


-1 








0-1 


) 


7i 







1 










V 


-1 / 








V 1 











/ 


(° 


i \ 








t 





-1 





\ 





-i 

















1 







-z 


) 


73 




1 













V * 


/ 








v 


-1 








/ 






/ 








-1 























-1 








75 


:= ^70717273 = 




-1 





















V 


- 


-1 





0^ 









Recall that: 

1. 7q = 7I = E4, the identity matrix; 

2. 7fe = -^4, fc = l,2,3; 

3- 7j7fe = ~lklj, j, k = 0, 1, 2, 3, 5, j ^ k. 

The transforms A and A" 1 have the following algebraic properties |5] 

1. A7i72737i-4~ 1 [p] = hp; 

2. A7i727372^~ 1 [p] = «2/o; 

3. A71 727373^^ [p] = -hp; 

4. A7i72737o-4 _1 [p] = (riii 

5. A7i727 3 A _1 [p] = -ipi 2 . 

The Dirac operator for a free massive particle of spin 1 /2 is 
D[$] := ^7o9 f - + [*]> 
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where m £ K and $ is a C 4 -valued function. Consider the quaternionic differ- 
ential operator 

R := Pf(id t + D) + P{ (-id t + D) - mM l \ 

where := \M^ 1±l1 ^ and M x [f] := /A, A e H(C). Using the algebraic 
properties of A and A^ 1 we obtain the following equality 

R = -A7i727 3 IIM~ 1 . (12) 

We consider a time-harmonic null-solution $ of the operator ID which has 
the form &(t,x) — q(x)e luJt , where u> <E M. and q(x) is a C 4 -valued function 
depending only on x. We have the following equation for q 

^, m q{x) = 0, 

where D w , m := ilo^q — Ylk=i Ikdk + iml. The corresponding quaternionic re- 
formulation of Ef w . m is the operator 

D a :=D + M a , 

where a := —{iujii + miz). That is 

D a = -A lll2 l^, m A- 1 . (13) 



3 Radiation condition for the Dirac operator 

In order to prove the Cauchy integral formula for null-solutions of the operator 
D u for the exterior domain K 3 \fi we need some appropriate radiation condition 
at infinity. 

Let us consider the equation 

D u K(x) = S(x), v^O, iel 3 . 

When Im v = it admits two solutions decreasing at infinity obtained (according 
to (Q)) by the application of the operator —D_ v to the fundamental solutions 
of the Helmholtz operator m ± (x) :— — e i7r \J respectively: 



/C ± (a;) = \ v + T iv T~\J ' u± ( x )- 

In order to omit one of these possibilities we impose the following radiation 
condition 

\v w + iu-rr I • fC(x) = o ( r— 7 ] , when \x\ — > oo. (14) 

I \x\ 2 \x\) \\x\J 
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It is easy to see that 



[ v 77 + iv-. — r | JC + (x) = O \ — 77 | , when \x\ — > oo 

V N Ny vw V 

and hence /C + fulfills (|l4|). For the function Kr this is not true. Note that K, + 
is precisely the fundamental solution JC V . 

We are now ready to prove the following 

Theorem 4 (Quaternionic Cauchy integral formula for the exterior domain) 
Let 

fe C* 1 (M 3 \n;H(C))nC(]R 3 \r2;H(C)), / e kerA,(K 3 \Tt), lmv>0 
and let f satisfy the radiation condition 

V m My Vfi/ 

XTien 

f(x) = -K u [f](x), Vxem. 3 \n. (15) 

Proof. Let T R be a sphere with center at the origin and radius R sufficiently 
large so that f2 is contained in the ball B R with boundary dB R = T R (Fig. ||). 

According to Theorem || in each point x of the domain il R := B R \ Q, we 
have the equality 



f{x) = / JC u (x-y)Tt(y)f(y)dT v 



y 
r 



/ K v (x-y)-^f(y)dT R 
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We now consider the limit of this equality when R — > oo. We have the following 
asymptotic relation 



Jr* \y\ Jr* \ \y\ \v\ / 



f(y)dT* i?^co. 



Using the radiation condition we obtain that this integral tends to zero when 
R —>■ oo. Thus, 

f(x) = J K v (x - y)Jt(y)f(y)dT y . 

■ 

The introduction of the radiation condition in the form ( |l4| ) allowed us 
to obtain a very simple proof of the Cauchy integral formula for the exterior 
domain. Nevertheless, as was shown in this condition can be represented in 
a more elegant form. We note that the expression which appears in (|l4|) can be 
rewritten as follows 

v[1 + £)+o(tK), W->oo. 




\x\; \\x\ 

Thus, a natural idea is to introduce the radiation condition as 



^ . + h) /w = °(h)' (16) 

because the term x/ \x\ apparently gives a faster decay. The problem here is 
that ( |l6| ) does not imply the decay of the product r% ■ f(x) when \x\ — > oo due 

to the fact that (1 + S) is a zero divisor. The following example makes this 

clear. Consider f{x) := (1 — S ) ■ \x\ 2 . This function obviously satisfies ( |l6| ) but 



f(x) = 0(\x\), \x\ — > oo. Thus, in order to be able to apply (16) instead of 



(pTD, we should prove that if / belongs to ker D l ,(R s \ tt) and satisfies (|16| ) then 
it decreases at infinity. Following arguments from |?j we prove this. 

Theorem 5 Let f € C 1 ^ 3 \ T!;H(C)) n C*(R 3 \ fi; 11(C)), / € ker A,(R 3 \ 
fl), Imp > 0, such that satisfies \l<\). Then 

[ \f(x)\ 2 c dr* = 0(l) as i?->oo. (17) 

J\x\=R 



Proof. Consider the expression 

2 



IX . 



= 2(|/|^-ImSc(/ — /) 
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Consequently, 



\x\=R 



fix) 



dT* = 2{ \f(x)\ldT* 

: J\x\=R 



(18) 



ImSc(/ f(x) — f{x)dT*)). 



x\=R 



Due to Theorem we have (using the notation of Theorem 0) 



/ { x ) — f{x)dT« = / / (x)Jt(x)f(x)dT x + 
x\=r \ x \ Jr 



n" 



(D r (f (x))-f(x) + f (x) ■ Df(x))dx. 



The first integral on the right-hand side is some constant C. In order to simplify 
the second we use Df = — vf and observe that 



D r f = -Df* = v*!* = v*f . 



Then 



(D r (f (x)) ■ f(x) + f (a;) • Df(x))dx = / iv*f (x) ■ f(x) ~ vf (x) ■ f(x))dx 
n R Jn R 



-2i1mv I f (x) ■ f(x)dx, 



we obtain that 
Sc 



/ ( x ) — f{x)dY" x \ =C -2ilmu / \f(x)\*dx 



\x\=R. 



Substituting this expression in ( pL8[ ) and using the radiation condition (16) we 
have 

/ \f(x)\ 2 c dr* = lmC -21mv [ \f(x)\ 2 c dx, R^oo. 
J\x\=r Jn R 

The last term above vanishes when Im^ = and is negative when Im^ > 0. In 
both cases ( p"7| ) is proved. ■ 

This theorem establishes the equivalence of the radiation conditions ( |l4] ) and 
( ^6| ) for the functions from ker D v . 

We now obtain the radiation condition for the Dirac operator in quaternionic 
form, that is for D a , and prove the Cauchy integral formula for the exterior 
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domain for this operator. After that, using the transform A we will rewrite 
these results for the operator D w , m . 

Let a be a purely vectorial complex quaternion a = ~ct, a (fc 6, ~ct 2 ^ and 
consider v := \/~cf 2 G C where the sign of the root is chosen in a way such that 
Im v > 0. 

Let 

K^[f]{x) := P+K4f](x) + p-K-4f](x) (19) 

where the operators K± u are denned by jlO| ) and P ± := ^M^ u± ^ are two 
mutually complementary projection operators in the set of H(C)-valued func- 
tions. 

In the case when a := ~ot G 6, v = 0, wc introduce the operator as 
follows 

Knf := K f - M^V f (20) 

where 

V [f](x) := /" o (z - y)Tt(y)f(y)dT y , x G K 3 \ r. (21) 
r 

Here # (£ —y) = — t- r — r- 

For all a € H(C) (that is a e 6 or a ^ 6), K-&f admits the following 
representation 

K^[f](x) = - [ {9 v (x - y) (-^^ iv ^-v) \ Tt(y)f(y) + 

J Wx-vv \ x -y\ ) 



+ e iJ {x-y)yt(y)f{y)-$}dT yi x G F. 



Theorem 6 Let f G C\R 3 \ fi; H(C)) n C(JR 3 \ Q; H(C)), / G kerD a (M 3 \ fi), 
a G H(C), Sea = such that f satisfies the radiation condition 

ix _ ( 1 \ 
vf {x) + -j — rf(x)~ct = o -. — r I , when \x\ — * oo, (22) 



X \ X 



where v := \J~ct 2 G C and Imis > 0. If a £ & then we suppose additionally that 
f(x) = o(l). Then 

f(x) = K- Bt [f\(x), VxeR 3 \Tl. 
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Proof. First we suppose that "57 &. Multiplying ( p2] ) by ~o? jv from the 
right-hand side we obtain 

m^+^m=o('\ (23) 



adding and subtracting (p2j) and (E3J) we obtain 



X± R (^Sl/W^ 



= (l±i)/(»)(,±^) = o( ] l). (24) 
Which can be writen as follows 

p+ (0 + r) H + p (( l - r) H = (r 1 when N - °°- 



(25) 



Thus, (22) is equivalent to fl2q), from which it can be seen that P + f fulfills the 
radiation condition ( |l6| ) and P~ f fulfills its conjugate which corresponds to the 
operator £>_„. Consequently the integrals K± u P ± f taken over the sphere T R 
(see the proof of Theorem ||) decrease when R — > oo. 

Since K± L , commute with P^ we obtain that the integral K-&f(x) taken 
over the sphere T R also decreases when R — ► oo. 

In the case when a is a zero divisor the radiation condition (E^) becomes 
X -f~c? = o ( - — r ) , when \x\ — > oo. (26) 



x \ x 



Since the behavior of 6q in ( pl| ) is of the type O when |a;| — > oo and 

since the expression M a Vq contains the multiplication from the right-hand side 

by ~at , and in this case fat = o (j^) > it can be seen that the integral ( |20| ) taken 

over the sphere T R decreases at infinity also. ■ 

We now consider the Dirac operator it := —{iuiix + mii). Note that 
~ct G S iff lo 2 = m 2 . In the following lines we will obtain the radiation condition 
for the Dirac operator in its classical representation using 7-matrices. 

Let us denote 

q := A" 1 /- (27) 
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According to (|l3|), if / e ker then q £ kerB w , m . Suppose that / fulfills the 
radiation condition ( |22| ) then we have 



Let us denote 
and consider the product 

xF = (xih + x 2 t2 + x 3 i 3 )F. 
Using properties 1.-3. of the transforms A and A' 1 (Section 2.3), we obtain 



xF = A717273 (^2 x k"/kj A 1 [F]. 



Consequently, using the properties 4. and 5., we have 



xF(x) = A717273 (^Y^Xk^fkj A 1 [f{x) ■ (-iu>h - mi 2 )] = 
= ^717273 (^2^k"fkj A^ 1 A-f 1727370 A^ 1 f(x) ■ (-ko)+ 



+A717273 ^Xk/fk A 1 yl7i727 3 A 1 f{x) ■ (-im) = 



\k=l 



Using (^7j) we can rewrite the previous equation as 

xF = A (^Y^Xklkj (-«^7o - im)q. 



Thus ( p2| ) can be rewritten as follows 

•'I 



i 1 ^ / 1 

"/ + TZ\f~& = vA 1 + A TZ\H x klk{u-lo + m)q = o(j— 



k=l 
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and hence the radiation condition for the Dirac operator in its traditional form 
is obtained in the following form 



vq(x) - (w7o - m)j^q(x) = o (^j > \ x \ ~> 00 ( 28 ) 

where := J2k=i x kJk- 

We introduce the following operator 

where ~ct = —{10111+17112). Using this operator, in pj an analogue of the Cauchy 
integral formula for the Dirac operator ID W)m in bounded domains was proved. 
Here with the aid of the obtained radiation condition we obtain this fact for the 
unbounded domains. 

Theorem 7 ( Cauchy integral formula for B Wim in the exterior domain R 3 \f2 y ) 
Let q £ C 1 (R 3 \fi;H(C)) n C(IR 3 \fi; H(C)), q £ ker© Wim (M 3 \fi) ; and q satisfies 
the radiation condition (\2ai). Then 



q(x) = K^ m [q](x), x£(R 3 \n). 



The proof consists in the application of (|13|) and Theorem ^. 
Acknowledgment: This work was partially supported by CONACYT, project 
32424-E. 
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